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Abstract 

By numerically solving the relativistic Boltzmann equations, we compute the time scale for relaxation to thermal equilibrium 
for an optically thick electron-positron plasma with baryon loading. We focus on the time scales of electromagnetic interactions. 
The coUisional integrals are obtained directly from the corresponding QED matrix elements. Thermalization time scales are 
computed for a wide range of values of both the total energy density (over 10 orders of magnitude) and of the baryonic 
loading parameter (over 6 orders of magnitude). This also allows us to study such interesting limiting cases as the almost 
purely electron-positron plasma or electron-proton plasma as well as intermediate cases. These results appear to be important 
both for laboratory experiments aimed at generating optically thick pair plasmas as well as for astrophysical models in which 
electron-positron pair plasmas play a relevant role. 

PACS numbers: 52.27.Ep,52.27.Ny,31.70.Hq 
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Current interest in electron-positron plasmas is due to 
the exciting possibility of generating such plasmas in lab- 
oratory facilities already operating or under construction, 
see e.g., [H 01, for a review see l3|]. Impressive progress 
made with ultra-intense lasers [J| has led to the creation 
of positrons at an unprecedented density of 10^^ cm""^ 
using ultra-intense short laser pulses, in a region of space 
with dimensions on the order of the Debye length. How- 
ever, such densities have not yet reached those neces- 
sary for the creation of an optically thick pair plasma 
@, 0] • Particle pairs are created at the focal point of 
ultra-intense lasers via the Bethe-Heitler conversion of 
hard x-ray bremsstrahlung photons [l] in the collision- 
less regime 0] . The approach to an optically thick phase 
may well be envisaged in the near future. 

Electron-positron plasmas are known to be present in 
compact astrophysical objects, leaving their characteris- 
tic imprint in the observed radiation spectra Op- 
tically thick electron-positron plasmas do indeed pjay a 
crucial role in the gamma-ray burst phenomenon [3|, |9| . 

From the theoretical point of view electron-positron 
pair plasmas are interesting because of the mass symme- 
try between the plasma components. This symmetry re- 
sults in the absence of both acoustic modes and Faraday 
rotation; waves and instabilities in such plasmas differ 
significantly from asymmetric electron-ion plasmas, see 
e.g. Besides, theoretical progress in understand- 

ing quark-gluon plasma in the high-temperature limit is 
linked to understanding QED plasma since the results in 
these two cases differ only by trivial factors containing 
the QCD degrees of freedom (color and flavor) . 

Most theoretical considerations so far have assumed 
that an electron-positron plasma is formed either in ther- 
mal equilibrium (common temperature, zero chemical po- 



tentials) or in chemical equilibrium (nonzero chemical 
potentials), see e.g. Q and references therein. However, 
it is necessary to establish the time scale for actually 
reaching such a configuration. The only way for parti- 
cles to thermalize, i.e., reach equilibrium distributions 
(Bose-Einstein or Fermi-Dirac) is via collisions. Colli- 
sions become relevant when the mean free path of the 
particles becomes smaller than the spatial dimensions of 
the plasma, and so the optical thickness condition is cru- 
cial for thermalization to occur. 

Thermalization (chemical equilibration) time scales for 
optically thick plasmas are estimated in the literature 
by order of magnitude arguments using essentially just 
the reaction rates of the dominant particle interaction 
processes, see e.g. ll|, |l2|. They have been computed 
using various approximations. In particular, electrons 
have been considered ultrarelativistic, and Coulomb log- 
arithm has been replaced by a constant. The accurate 
determination of such time scales as presented here is in- 
stead accomplished by solving the relativistic Boltzmann 
equations including the collisional integrals representing 
all possible particle interactions. In this case the Boltz- 
mann equations become highly nonlinear coupled partial 
integro-differential equations which can only be solved 
numerically. 

We developed a relativistic kinetic code treating the 
plasma as homogeneous and isotropic and have previ- 
ously determined the thermalization time scales for an 
electron-positron plasma for selected initial conditions 



13 1 . This approach was generalized to include protons in 
14| . We focus only on the electromagnetic interactions, 
which have a time scale of less than 10 sec for our sys- 
tem, and therefore on the proton and leptonic component 
of the plasma. The presence of neutrons and their possi- 
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ble equilibrium due to weak interactions will occur only 
on much longer time scales. 

In this paper we report on the systematic results ob- 
tained by exploring the large parameter space character- 
izing pair plasmas with baryonic loading. The two basic 
parameters are the total energy density p and the bary- 
onic loading parameter 
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where ph and pe,7 are respectively the total energy densi- 
ties of baryons and electron-positron-photon plasma, Up 
and Trip are the proton number density and proton mass, 
and c is the speed of light. We choose the following range 
of plasma parameters 



10^^ < p< 10^^ erg/cm^, 
10"^ <B< 10^, 



(2) 
(3) 



allowing us to also treat the limiting cases of almost pure 
electron-positron plasma with i? <iC 1, and almost pure 
electron-ion plasma with B ~ mp/rrie, respectively. The 
temperatures in thermal equilibrium corresponding to ([2]) 
are 0.1 < < 10 MeV. 

Given the smallness of the plasma parameter g = 
(rieAp)^^ <^ 1, where Aq is the Debye length and 
is the electron number density, it is sufficient to use 
one-particle distribution functions. In fact, for the pure 
electron-positron plasma, the inequality 3 10^^ < 5 < 
10~^ holds in the region of the temperatures of interest. 
In a homogeneous and isotropic plasma the distribution 
functions /(e,i) depend on the energy e of the particle 
and on the time t. We treat the plasma as nondegenerate, 
neglecting neutrino channels as well as the creation and 
annihilation of baryons and the weak interactions 14 1 . 

The relativistic Boltzmann equations 15|, |16| for pho- 
tons, electrons, positrons, and protons in our case are 



c dt 



(4) 



where the index i denotes the type of particle and jyf, xf 
are the emission and the absorption coefficients for the 
production of the ith-particle via the reaction labeled by 
q. We account for all relevant binary and triple interac- 
tions between electrons, positrons, photons, and protons 
as summarized in Tables HI and HIl 

It has been shown [l^] that independent of the func- 
tional form of the initial distribution functions fi{e,0), 
plasma evolves to a thermal equilibrium state through 
the kinetic equilibrium, when the distribution functions 
of all the particles acquire the same form 



/,(e) = cxp - 



Binary interactions 


Radiative and 
pair producing variants 


M0ller and Bhabha 

^1 ^2 ^ ^1 ^2 


Bremsstrahlung 


Single Compton 

e*7 — ye^-y' 


Double Compton 
e='=7-f->e*'7'7" 


Pair production 
and annihilation 


Radiative pair production 
and 3-photon annihilation 

77'-H-e='=e^7" 
e''^e^-f^77'7" 




e±7^e±'e=Fe±" 



TABLE I: Microphysical processes in the pair plasma. 



Binary interactions 
(Coulomb scattering) 


Radiative and 
pair producing variants 


PiP2 — y p'iP'2 


± ' ±1 
pe -H- p e 7 


pe — > p e 


pj-^r^p'e^e^ 



TABLE II: Microphysical processes in the pair plasma involv- 
ing protons. For details see also Q|. 



where = ei/{miC^) is the energy of the particles, 
Lpi = pi/{miC^) and 6i = k-sTi / {rriiC^) are their chemical 
potentials and temperatures, and fee is Boltzmann's con- 
stant. The unique signature of kinetic equilibrium is the 
equal temperatures of all the particles and the nonzero 
chemical potential of the photons. In fact the same is 
also true for a pair plasma with proton loading [l3| . The 
approach to complete thermal equilibrium is more com- 
plicated in this latter case and depends on the baryon 
loading. For B ^ y^nip/me, protons are rare and ther- 
malize via proton-electron (positron) elastic scattering, 
while in the opposite case B ^ ^ mp/rrif,, proton-proton 
Coulomb scattering dominates over the proton-electron 
scattering and brings protons into thermal equilibrium 
first with themselves. Then protons thermalize with the 
pair plasma through triple interactions, for details see 
[ij] . The two-body time scales involving protons should 
be compared with the three-body time scales bringing 
the electron-positron-photon plasma into thermal equi- 
librium. In fact we found that for S <C 1, the electron- 
positron-photon plasma reaches thermal equilibrium at 
a given temperature, while protons reach thermal equi- 
librium with themselves at a different temperature; only 
later the plasma evolves to complete thermal equilibrium 
with the single temperature on a time scale 



(5) 



Tth ^ Max [rsp, Min [r^p, Tpp)] 



(6) 
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FIG. 1: The thermalization time scale of the electron- 
positron-photon component of plasma as a function of the 
total energy density and the baryonic loading parameter. The 
energy density is measured in erg/cm"^, time is seconds. 



FIG. 2: The final thermalization time scale of a pair plasma 
with baryonic loading as a function of the total energy density 
and the baryonic loading parameter. The energy density is 
measured in erg/cm"^, time is seconds. 



where 




I lip I ^-l 

me 

T3p ~ {aaTUecy^ 



(7) 

(8) 
(9) 



are the proton-electron (positron) clastic scattering time 
scale, the proton-proton elastic scattering time scale, 
and the three-particle interaction time scale respectively, 
while ctt is the Thomson cross-section and a is the fine 
structure constant. In Q-lIHl) the energy dependence of 
the corresponding time scales is neglected. 

The chemical relaxation (thermalization) time scale is 
usually computed as 



(10) 



where Fi — exp{(pi/6i) is the fugacity of a particle of 
type i. Instead of Fi we use one of the quantities 9i, ipi, 
rii, or Pi in this computation. 

We solved the Bohzmann equations with parameters 
{p,B) in the range given by Eqs. ([2]) and dS]). In total 
78 models were computed, starting from a nonequilib- 
rium configuration until reaching a steady state solution 
on the computational grid with 20 intervals for the par- 
ticle energy and 16 intervals for the angles, for details 
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FIG. 3: The final thermalization time scale of pair plasma 
with baryonic loading as a function of the total energy den- 
sity for selected values of the baryonic loading parameter 
B = (10"^, 10"^•^ 1, 10, 10^ 10^). The energy density is mea- 
sured in erg/cm'^, time is seconds. Error bars correspond to 
one standard deviation of the time scale (|lll) away from the 
average value Tth over the interval tin < t < tfi„. 



see [14| . For each model we computed the corresponding 
time scales for all particles of the zth kind. For practical 
purposes, instead of (|10|) we used the following approxi- 
mation 



1 



Tth 



fin 



''"[0(t)-0(W)](^) dt, (11) 



with tin < if in < imax, whcre tmax is the moment of time 
where the steady solution is reached and tin and t fin are 
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the boundaries of the time interval over which the averag- 
ing is performed, for details see 17|. The thermalization 
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FIG. 4: The thermalization time scale of the electron- 
positron-photon component of the plasma as a function of 
the total energy density (points), compared with the rgp time 
scale (joined points) computed using © for B = 1. The 
energy density is measured in erg/cm'^, time is seconds. 
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FIG. 5: The final thermalization time scale of a pair plasma 
with baryonic loading as a function of the total energy den- 
sity (points), compared with the Tth time scale (joined points) 
computed using (|6} for B = 1. The energy density is mea- 
sured in erg/cm'', time is seconds. 



time scale of the electron-positron-photon component is 
shown in Fig. [l]as a function of the total energy density 
of the plasma and the baryonic loading parameter. The 
time scales of electrons, positrons and photons coincide. 
The final thermalization time scale of pair plasma with 
baryonic loading is shown in Fig. [5] Its dependence on 
either variable cannot be fit by a simple power law, al- 
though it decreases monotonically with increasing total 
energy density, while it is not even a monotonic function 
of the baryonic loading parameter. 

In Fig.[3]the final thermalization time scale is shown for 
all the models we computed, along with the "error bars" 
which mark one standard deviation of the time scale pT|) 
away from the average value Tth hi the averaging interval 
tin < t < tfin- The largest source of error comes from 
the small values of the time derivative in dTll . although 
errors are typically below a few percent. 



In Fig. |4] we compare for B = 1 the actual value of the 
thermalization time scale of the electron-positron-photon 
component with the value estimated from ([9]). Both val- 
ues clearly differ significantly. Actually the systematic 
underestimation by more than one order of magnitude 
which occurs for _B < 1 disappears for larger baryonic 
loading. 

In Fig. [5] we present the computed values of the final 
thermalization time scale of the pair plasma with bary- 
onic loading together with the value estimated from ([6|), 
again for B = 1. Unlike the previous case, the final ther- 
malization time scale is a more complex function of the 
total energy density. Interestingly, less significant devi- 
ations from the value (|6]) occur at the extremes of the 
interval (|3]). 

In this paper we have computed for the first time 
the time scale of thermalization for an electron-positron 
plasma with proton loading over wide ranges of both 
the total energy density (10 orders of magnitude) and 
baryonic loading parameter (6 orders of magnitude) al- 
lowing the treatment of the limiting cases of almost 
pure electron-positron plasma, almost pure electron-ion 
plasma as well as intermediate cases. The final result 
is presented in Fig. [T] and [2l The relaxation to thermal 
equilibrium for the total energy density ([2|) always occurs 
on a time scale less than 1Q~^ sec. It is interesting that 
the electron-positron-photon component and/or proton 
component can thermalize earlier than the time at which 
complete thermal equilibrium is reached. The relevant 
time scales are given and compared with the ordcr-of- 
magnitude estimates. Unlike previous work there are no 
simplifying assumptions in our method since collisional 
integrals in the Boltzmann equations are computed di- 
rectly from the corresponding QED matrix elements, e.g. 
from the first principles. 

These results may be of relevance for the ongoing and 
future laboratory experiments aimed at creating electron- 
positron plasmas. Current optical lasers producing pulses 
during ^ 10~^^ sec carrying energy ~ 10^ J= 10^ erg 
are capable to produce positrons with the number den- 
sity 10^^ cm^"^ 0. There are claims that densities of 
the order of 10^^ cm~^ are reachable [l3|- These den- 
sities today are yet far from 10^^ cm~^ required for the 
plasma with the size ro — nm to be optically thick [Bj. 
Notice, that the expansion timescale of such plasma will 
be tq/c ^ 10~^* sec, while the timescale to establish ki- 
netic equilibrium for the number density considered is 
of the same order of magnitude. These arguments show 
that theoretical results obtained assuming thermal or ki- 
netic equilibrium, such as in Q, cannot be applied to 
pair plasma, generated by ultraintense lasers. 

However, results presented in this paper are impor- 
tant for understanding astrophysical systems observed 
today in which optically thick electron-positron plasmas 
are present. As specific example we recall that electron- 
positron pairs play the crucial rule in the dynamics of 
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GRB sources. Considerin g ty pical energies and initial 
radii for GRB progenitors [19| 

lO^'^erg <Eo < lO^^erg, let'em < Rq < lO^cm, (12) 

we estimate the range for the energy density in GRB 
sources 



10^3^ <p< 10^2^, 
cm^ cm'^ 



(13) 



which coincides with ([2]). As for the baryonic loading of 
GRBs it is typically in the lower range of ([2|) , namely 



10"^ < B < IQ-^. 



(14) 



Such high energy density leads to large number density 
of electron-positron pairs in the source of GRB, of the 
order of 



10^° cm"^ <n< 10^^ cm"^, 



(15) 



making it opaque to photons with huge optical depth of 
the order of 



10" <T < 10". 



(16) 



In fact, the radiative pressure of optically thick electron- 
positron plasma in these systems is resp onsible for the 
effect of accelerated expansion [20l422l |. leading to un- 



precedented Lorentz factors attained F 
10'^, see e.g. 



B ^, up to 

23l , |24| . The role of the baryon admixture 



in electron-positron plasma in GRBs is to transfer inter- 
nal energy of pairs and photons into kinetic energy of 
the bulk motion thus giving origin to afterglows of GRBs 
[l^ . Notice that in GRBs the timescales of thermal- 
ization are much shorter than the dynamical timescales 
Rq/c ~ 10"'^ sec, which implies that expanding electron- 
positron plasma even in the presence of baryons is in ther- 
mal equilibrium during the accelerating optically thick 
phase [255. 

After completion of this work we learned about the 
publication of [2^ where work similar to ours has been 
performed. Between this paper and our work conceptual 
differences should be noted which concern the attribution 
of thermalization to two-body M0ller and Bhabha scat- 
tering, while we have pointed out explicitly that three- 
body interactions play an essential role. The thermaliza- 
tion time scales obtained by us have been computed with 
reference to these three-body interactions. 



[1] J. Myatt et al., Phys. Rev. E 79, 066409 (2009). 
[2] M. H. Thoma, Reviews of Modern Physics 81, 959 
(2009). 

[3] R. RufEni, G. Vereshchagin, and S.-S. Xue, 

Phys. Rep. 487, 1 (2010). 
[4] H. Chen et al., Phys. Rev. Lett. 102, 105001 (2009). 
[5] J. I. Katz, ApJS 127, 371 (2000). 

[6] M. G. Mustafa and B. Kampfer, Phys. Rev. A 79, 020103 
(2009). 

[7] S. C. Wilks, W. L. Kruer, M. Tabak, and A. B. Langdon, 
Physical Review Letters 69, 1383 (1992). 

[8] E. Churazov, R. Sunyaev, S. Sazonov, M. Revnivtsev, 
and D. Varshalovich, MNRAS 357, 1377 (2005). 

[9] R. RufEni et al.. The Blackholic energy and the canon- 
ical Gamma-Ray Burst IV: the "long," "genuine short" 
and "fake-disguised short" GRBs, in American Institute 
of Physics Conference Series, edited by M. Novello and 
S. Perez, volume 1132 of American Institute of Physics 
Conference Series, pages 199-266, 2009. 
[10] G. P. Zank and R. G. Greaves, Phys. Rev. E 51, 6079 
(1995). 

[11] R. J. Gould, Physics of Fluids 24, 102 (1981). 

[12] S. Stepney, MNRAS 202, 467 (1983). 

[13] A. G. Aksenov, R. Rufiini, and G. V. Vereshchagin, 
Phys. Rev. Lett. 99, 125003 (2007). 

[14] A. G. Aksenov, R. Ruffini, and G. V. Vereshchagin, 
Phys. Rev. D 79, 043008 (2009). 

[15] S. Belyaev and G. Budker, DAN SSSR 107, 807 (1956). 

[16] D. Mihalas and B. W. Mihalas, Foundations of Radia- 
tion Hydrodynamics, New York, Oxford University Press, 
1984. 

[17] A. G. Aksenov, R. Ruffini, and G. V. Vereshchagin, AIP 

Conf. Proc. 1205, 11 (2010). 
[18] B. Shen and J. Meyer-Ter-Vehn, Phys. Rev. E 65, 016405 

(2002). 

[19] T. Piran, Phys. Rep. 314, 575 (1999). 

[20] R. Ruffini, J. D. Salmonson, J. R. Wilson, and S.-S. Xue, 

A&A 350, 334 (1999). 
[21] R. Ruffini, J. D. Salmonson, J. R. Wilson, and S.-S. Xue, 

A&A 359, 855 (2000). 
[22] C. L. Bianco, R. Ruffini, G. Vereshchagin, and S.-S. Xue, 

J.Korean Phys.Soc. 49, 722 (2006). 
[23] A. A. Abdo et al.. Science 323, 1688 (2009). 
[24] L. Izzo et al., in Proceedings of the First Galileo - Xu 

Cuangqi meeting held in October 26-30, 2009 Shanghai 

(China), in press, 2010. 
[25] A. G. Aksenov, C. L. Bianco, R. Ruffini, and G. V. 

Vereshchagin, AIP Conf. Proc. 1000, 309 (2008). 
[26] I. Kuznetsova, D. Habs, and J. Rafelski, ArXiv e-print 

arXiv0911.0118 (2009). 



Acknowledgements. We thank both anonymous ref- 
erees for their comments which allowed to improve re- 
markably the paper. 



5 



